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Abstract 

Accurately calibrated effective field theories are used to compute atomic parity non-conserving 
(APNC) observables. Although accurately calibrated, these effective field theories predict a large 
spread in the neutron skin of heavy nuclei. While the neutron skin is strongly correlated to a 
large number of physical observables, in this contribution we focus on its impact on new physics 
through APNC observables. The addition of an isoscalar-isovector coupling constant to the ef- 
fective Lagrangian generates a wide range of values for the neutron skin of heavy nuclei without 
compromising the success of the model in reproducing well constrained nuclear observables. Earlier 
studies have suggested that the use of isotopic ratios of APNC observables may eliminate their 
sensitivity to atomic structure. This leaves nuclear structure uncertainties as the main impediment 
for identifying physics beyond the standard model. We establish that uncertainties in the neutron 
skin of heavy nuclei are at present too large to measure isotopic ratios to better than the 0.1% 
accuracy required to test the standard model. However, we argue that such uncertainties will be 
significantly reduced by the upcoming measurement of the neutron radius in 208 Pb at the Jefferson 
Laboratory. 

PACS numbers: 24.80. +y, 21.10.Gv, 32.80.Ys, 11.30.Er 
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I. INTRODUCTION 



The quest for discovering fundamental characteristics of atomic nuclei, such as their sizes 
and shapes, is as old as Nuclear Physics itself. Although for the most part this challenge has 
been met with enormous success, precise knowledge of certain fundamental nuclear-structure 
properties — among them the spatial distribution of neutrons — is still lacking. While the 
original motivation for measuring the neutron density of atomic nuclei was deeply rooted in 
nuclear structure, the quest has been recently revived in response to the widespread impact 
that such a measurement will have over seemingly unrelated areas of physics. Even though 
the present study concentrates exclusively on atomic parity non-conservation (APNC) 0, 0| , 
a brief review of areas that could benefit from an accurate measurement of the neutron 
density is also presented. 

Unquestionably, the accurate determination of the neutron density — and in particular its 
root-mean-square (rms) radius — remains a top priority in nuclear-structure physics. This in- 
adequate state of affairs is in stark contrast to our present knowledge of proton densities. In- 
deed, the proton distribution in nuclei has already been mapped with exquisite precision, of- 
ten much better than 1%! , due to the availability of state-of-the-art electron-scattering fa- 
cilities all over the world. A particularly illustrative example is 208 Pb, a nucleus whose charge 
radius is presently known with an extremely high accuracy (r c h = 5.5010 ± 0.0009 fm |3|) 
that exceeds at least by two orders of magnitude that of its neutron radius. Clearly, an ac- 
curate determination of neutron radii is a pressing issue. Further, without such knowledge 
for stable nuclei, the prospects of real progress in the future domain of nuclear structure, 
namely, exotic nuclei, may be seriously compromised. 

Precision studies of the nucleon-nucleon (NN) force is an area that has been directly 
affected by our poor knowledge of neutron radii. At the most fundamental level, Quantum 
Chromodynamics predicts a small flavor violation in the strong force due to the different 
charges and masses of the up and down quarks. In turn, this flavor violation induces a 
small (of the order of 1%) breaking in the isospin symmetry of the NN force. However, to 
accurately quantify these "novel" isospin violations, it is necessary to have all "conventional" 
effects under control. A good example of this interplay between novel and conventional is 
provided by the Nolen-Schiffer (or Coulomb energy) anomaly. The anomaly consists in the 
residual discrepancy observed in the ground-state energies of mirror nuclei with a single 
nucleon added (or removed) from a closed shell (e.g., 41 Sc- 41 Ca). However, some of the 
conventional effects depend sensitively on our knowledge of neutron radii. Indeed, it has 
been argued that by using neutron skins — defined as the difference between the neutron 
and proton rms radii — significantly smaller than those predicted by theory (and apparently 
justified by experiment) the anomaly disappears 

Knowledge of neutron radii could also remarkably improve our understanding of the 
equation of state of neutron-rich matter. Although existing ground-state observables appear 
to constrain the symmetry energy of nuclear matter at a neutron density around p n = 
0.10 fm~ 3 (as a reference value typical in nuclei) |H E|, the density dependence of the 
symmetry energy is unknown. Recently, however, it has been found that the neutron skin of 
a heavy nucleus, like 208 Pb, calculated with different non-relativistic and relativistic mean 
field parametrizations displays a tight linear relationship with the slope of the equation 
of state of neutron matter (or, analogously, with the derivative of the symmetry energy), 
evaluated at p n ~ 0.10 fm~ 3 jll El 0, la El • in particular, models with a "stiffer" symmetry 
energy predict larger neutron skins. For instance, the neutron skin in 208 Pb computed with 
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relativistic models is larger by typically 0.1-0.2 fm than the value obtained with Skyrme 
forces. In turn, the neutron skin of 208 Pb is linearly correlated with the neutron skin of 
other relatively heavy nuclei such as 132 Sn or 138 Ba. Thus, a measurement of the neutron 
radius of 208 Pb, or indeed of any other heavy nucleus, could fix a fundamental property of 
the equation of state. 

Fixing the slope of the symmetry energy at p = 0.10 fm -3 will also impact favorably 
on astrophysical observables, particularly on the structure and dynamics of neutron stars. 
Indeed, strong correlations between the neutron skin of 208 Pb and various neutron-star 
properties have been established ^(j. Among these, a model-independent relation between 
the neutron skin of 208 Pb and the transition density from uniform neutron-rich matter (in the 
mantle) to a non-uniform phase (in the crust) was observed. Other neutron-star properties 
sensitive to the neutron skin of 208 Pb include the radius of the star, its composition, and its 
cooling mechanism (ill 12l | . Ultimately, these correlations emerge as a result of the similar 



composition of the neutron skin of a heavy nucleus and the crust of a neutron star, namely, 
neutron-rich matter at similar densities. 

While the above arguments are already compelling enough to justify the commission of 
new experiments to measure the neutron radius of 208 Pb with unprecedented accuracy [13], 
this manuscript will focus on the crucial role that such a measurement could have on pre- 
cision studies of the standard model, vis-a-vis, atomic parity non-conservation. The road 
to new physics beyond the standard model in low-energy tests in atoms, which stems from 
violations of fundamental symmetries that occur in the weak interaction, passes through the 
observation of deviations between the values measured in the laboratory and the predictions 
of the standard model. The effects are inherently small and, traditionally, the analysis of 
APNC experiments has been hindered by uncertainties in both atomic and nuclear-structure 
theory. However, a fruitful strategy has been devised to remove the sensitivity to atomic 
theory. This strategy consists in measuring ratios of parity-violating observables along an 
isotopic chain ^]|. In this manner, uncertainties in the atomic theory factor out from the 



ratios. Based on their long chains of several naturally occurring isotopes, cesium, barium, 
dysprosium, and ytterbium appear as ideal candidates 0, 0, Ell, [3, 17, [l8j]. As a result, 
nuclear-structure uncertainties — primarily in the form of neutron radii — remain as the lim- 
iting factor in the search for physics beyond the standard model jlil 

In an earlier study, a strong correlation between the neutron radius of 208 Pb and the neu- 
tron radii of 138 Ba, 158 Dy, and 176 Yb was established [23|. In that study the density depen- 
dence of the symmetry energy — and consequently the neutron radius of 208 Pb — was modified 
through the addition of a new coupling constant (A v ) to the underlying Lagrangian 



This new coupling allows one to modify the density dependence of the symmetry energy 
without changing the saturation properties of symmetric nuclear matter. In a heavy nucleus 
A v modifies the rms radius of the neutron density while leaving the rms radius of the proton 
density and the total binding energy practically unchanged |10|,lll(. However, the nuclear- 
structure model used lacked both deformation and pairing effects that are important for the 
nuclei considered in APNC experiments. Thus, in Ref. |23[ the calculations were limited to 
the study of a single member of each isotopic chain (the one with a closed neutron shell). 
This shortcoming is resolved in the present contribution. As a result, one is now able to map 
the neutron radii of these long isotopic chains as a function of both the neutron number N 
and A v . There are previous studies of the isotopic dependence of the neutron radius in the 
relativistic mean field (RMF) approximation 2M, HI and of its dependence on A v |23|. Yet 
the merit of the present work is that uncertainties in the neutron radius are studied for the 
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first time as a function of both N and A v within a unified "best-fit" model. In addition, it 
has been recently pointed out that the accuracy in the measurements of APNC effects can 
be improved by using high-Z atoms, such as francium, where the parity violating effects are 
expected to be an order of magnitude larger than in cesium [26(. Thus we also include, in 
anticipation of future experiments, the analysis of the nuclear structure corrections to the 
weak charge in francium isotopes. 

The manuscript has been organized as follows. In Sec. |H] a short description of the 
nuclear model is given, with special emphasis paid to the treatment of deformation and 
pairing correlations as well as to the predictions of the model applied to cesium, barium, 
dysprosium, and ytterbium. In Sec. IIHI a brief review of atomic parity violation is given, 
with the main goal of addressing uncertainties in APNC observables emerging from our 
poor knowledge of neutron radii; a detailed analysis of such uncertainties within our model 
is also presented in the same section. Finally the conclusions are reserved for Sec. II VI We 
may summarize the two main conclusions of our work as follows. First, it appears that even 
when using the wealth of existing ground-state observables to constrain our nuclear-structure 
models, the uncertainties in the neutron radius of heavy nuclei — and consequently on APNC 
observables — may have been underestimated. Second, we have found, as others have before 
us, a tight correlation between the neutron radius of 208 Pb and the neutron radii of barium, 
dysprosium, and ytterbium, which is shown here to hold also in deformed systems and in 
the presence of pairing. This correlation suggests that the nuclear structure observables 
relevant for APNC will profit positively from the upco ming high precision measurement of 
the neutron radius of 208 Pb at the Jefferson Laboratory 13]. 

II. THE NUCLEAR MODEL 

The mean field treatment of effective field theories of hadrons, generally known as quan- 
tum hadrodynamics (QHD), is well established as a successful approach for describin g div erse 
bulk and single-particle properties of finite nuclei and uniform nuclear matter |27l . l28l . 29] . 
These effective field theories are based on a Lagrangian density that contains the nucleon (as 
an elementary Dirac particle) together with an isoscalar-scalar (a) meson, an isoscalar-vector 
(u) meson, an isovector-vector (p) meson, and the photon as the relevant degrees of freedom 
for describing the nuclear many-body problem. In the mean field approximation the me- 
son fields are replaced by their ground-state expectation values, thereby becoming classical 
fields. The quantum structure of the theory is carried by the nucleon field. For systems with 
time reversal symmetry, only the time-like component of the vector fields contribute. At the 
mean field ( "tadpole" ) level, charge conservation implies that only the third component (in 
isospin space) of the isovector p-meson field does not vanish. As a final product, one obtains 
the following Hamiltonian density: 
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where the summation runs over the occupied nucleon states <^ a (r) of positive energy and r 3 
stands for the third component of the isospin operator. The scaled meson fields associated 
with the a, u>, and p mesons are, respectively, $ = g s (po(r), W = g v Vo(r), and B = g p bo(r). 
Finally, A = eA Q (r) represents the time-like component of the photon field. 

The variation of Eq. (J1J with respect to the Dirac spinors yields the Dirac equation 
satisfied by the nucleons, while the variations with respect to the meson fields lead to the 
Klein-Gordon equations for the corresponding mesons [27ll28l l29|. The functional in Eq. 
contains cubic and quartic scalar meson self-interactions (couplings k and A, respectively) 
that are tuned to bring the value of the compression modulus of symmetric nuclear matter 
down to the empirical value (K = 200-300 MeV) |3(j. Further, through the inclusion of 
these terms one can accurately reproduce the systematics of finite nuclei throughout the 
periodic table as, e.g., in the case of the celebrated relativistic mean field (RMF) parameter 
set NL3 [3l|. The additional isoscalar-isovector cross coupling A v enables one to modify the 
density dependence of the symmetry energy without compromising the success of the model. 

Insofar as we are interested in the properties of some relatively heavy nuclei that have been 
identified as possible candidates for APNC experiments, we extend the model of Ref. [23| 
to include nuclear deformation and pairing correlations in the present RMF study. Only 
quadrupole deformation (/3 2 ) is included though, as it is known that hexadecapole deforma- 
tion is generally small for the nuclei under consideration 24]|. A BCS framework is employed 
to deal with the pairing correlations. Note that because the nuclei to be studied lie relatively 
close to the /3-stability line, the use of the BCS approximation is justified as the Fermi level 
is sufficiently bound and the mixing with continuum states — essential for drip line nuclei — is 
not relevant here. As in most BCS calculations, the strength of the neutron-neutron and 
proton-proton pairing interactions is determined from the neutron and proton pairing gaps 
(respectively, A n and A p ) evaluated from the experimental odd-even mass differences, using 
a five-point formula |32j. 

The Dirac equation and the Klein-Gordon equations for the meson fields are treated using 
an axially deformed basis expansion method as described in detail in Refs. [33l 0|- The 
calculations presented here have been carried out by expanding the fermionic wave functions 
and the bosonic fields in 12 and 20 oscillator shells, respectively. The solution of the equa- 
tions in an axially deformed basis for the odd-even and odd-odd nuclei is more complicated. 
Indeed, in these cases time reversal symmetry is broken and the odd particle induces polar- 
ization currents and time-odd components in the mean fields. However, the impact of these 
effects on deformation and binding energies is small [35| and will be neglected henceforth. In 
the pairing calculation of the odd nucleon we use the blocking approximation |36j], thereby 
restoring time reversal invar iance in the intrinsic frame. 

For the isotopic chains under consideration for APNC experiments, namely, Cs, Ba, Dy, 
and Yb, the NL3 parameter set yields ground-state properties, such as binding energies 
and deformations, in good agreement with the available experimental data j24[. For our 
calculations we adopt the NL3 parameter set and suitably modify it with the addition 
of the isoscalar-isovector coupling A v , prompted by the fact that the slope of the nuclear 
symmetry energy at (or near) the saturation density is at present unknown. Thus, we utilize 
the nonlinear isoscalar-isovector term A v to change the density dependence of the symmetry 
energy, which in turn modifies the neutron radius of neutron-rich nuclei. This can be done 
without a significant change in those ground-state properties of finite nuclei that are well 
constrained experimentally [lOj]. In practice, for a given A v value we readjust the nucleon- 
rho coupling constant g p so that the symmetry energy of nuclear matter at a density of 
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p = 0.10 far 3 be 25.68 MeV 0, Q 0. The choice is motivated by the fact that the 
symmetry energy of uniform nuclear matter at saturation density is not well constrained by 
the known properties of finite nuclei. Rather, it is some average between the bulk symmetry 
energy at saturation and the surface symmetry energy which is constrained by the binding 
energy of finite nuclei 0, 0|. For example, typical effective nuclear forces fitted to the 
empirical energies of nuclei predict a liquid drop model symmetry energy coefficient, which 
contains a surface correction, of 22-24 MeV for 208 Pb |37[ , far smaller than the value of the 
bulk symmetry energy coefficient for that interaction. 

Although it has been established that a variety of ground-state properties of spherical 



nuclei are insensitive to the choice of A v 1 10, 11, 231, this is not the case for the neutron 



skin. The neutron skin t = r n — r p , where r n and r p are the neutron and proton rms radii 
(r n (p) = ( r2 )n(p))) is strongly sensitive to the density dependence of the symmetry energy 
and hence, to the isoscalar-isovector coupling A v . This is because the slope (i.e., pressure) 
of the symmetry energy controls how far is the neutron skin pushed out relative to the 
symmetric core. In order to show that this scenario holds true even when the spherical 
symmetry is broken and pairing correlations are important, we display in Fig. ^ the charge 
radius, binding energy, quadrupole deformation, and the neutron skin as a function of A v 
for the representative heavy nucleus 174 Yb. Note that all observables in Fig. ^ have been 
normalized to the corresponding NL3 values (A v = 0). It is seen that the charge radius, 
binding energy and quadrupole deformation are insensitive to the value of A v in the 0-0.025 
range; the quadrupole deformation displays the largest change, but this amounts to only 2%. 
In contrast, the neutron skin decreases rapidly with increasing A v and it varies by about 
25% when A v is changed from to 0.025. This reflects the softening of the symmetry energy 
as A v increases. Therefore, we conclude that neither pairing correlations nor deformation 
alter our earlier conclusions about the role of A v in the systematics of finite nuclei. 

The neutron skin of a nucleus is an important component in the study of atomic parity 
non-conservation. So is its charge radius. While neutron radii of heavy nuclei (r n ) are 
poorly known, high precision data for the charge radii are available. To test the reliability 
of our model, we now compare in Fig. |2] calculated charge radii with experimental data 
for some of the relevant isotopes used in APNC experiments; these include Cs, Ba, Dy 
and Yb. The experimental values are taken from the recent compilation of Ref. [3j], which 
combines measured data from electron scattering, from muonic X-rays, and from K a and 
optical isotope shifts for the purpose of providing a unified set of nuclear charge radii. 
Theoretical calculations are presented with the original NL3 parametrization (A v = 0) and 
with one for which the density dependence of the symmetry energy has been softened by 
fixing the value of the isoscalar-isovector coupling to A v = 0.025. It can be seen that both 
sets of calculated charge radii are in excellent agreement with the experimental values for 
all of the considered nuclei. At a more quantitative level, we find that the relative x 2 value 
between the theoretical and experimental charge radii for the isotopic chains of Cs, Ba, Dy, 
and Yb is always smaller than 0.25% for both (A v = and A v = 0.025) parameter sets. 
This result validates the use of both parameter sets in computing key nuclear observables 
of relevance to the APNC program. Note that the binding energies and deformations for 
all the nuclei considered here are also in good agreement with the experimental data (not 
shown). Moreover, and as alluded earlier, these observables are highly insensitive to the 
choice of A v . 

A more daring challenge for a theoretical model than the charge radii themselves is the 
reproduction of the observed charge isotope shifts. One famous case is the pronounced kink 
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exhibited by the experimental charge radii of the Pb nuclei at the magic neutron number 
N = 126. The kink is not reproduced by conventional non-relativistic theories (either 
Skyrme or Gogny forces). In contrast, the RMF calculations show a remarkable agreement 
with experiment j28[. The origin for the deviations between the two models has been traced 
back, largely, to the isospin dependence of the spin-orbit term. We present in Fig. El results for 
the isotope shifts in two of our flagship chains for the subsequent APNC analysis: the cesium 
and francium alkali metals. The shifts are referred to 133 Cs and 212 Fr, respectively, as are 
the data of Ref . 0| ■ The experimental method of preference for the measurement of nuclear 
charge isotope shifts is high-resolution laser spectroscopy, which exploits the isotopic effects 
on the hyperfine structure of atomic optical transitions. In addition to the RMF calculations 
with A v = and A v = 0.025 we have drawn in Fig. El the SkM* and SIII results of the non- 
relativistic Hartree-Fock calculations of Ref. jij for Cs. One immediately recognizes the 
ability of the RMF model in reproducing the trend of the changes of the experimental data 
when neutrons are added; conspicuously, the change of slope at the N = 82 shell closure for 
Cs and at the N = 126 shell closure for Fr. In spite of this, the theoretical predictions may 
be off the experiment by a factor of ~ 2 in some specific isotopes. One should keep in mind 
that this happens in cases where the isotope shift is a small quantity and that it emerges 
from the cancellation of two much larger numbers (the square charge radii of two isotopes). 
Similarly to the observables discussed above, the calculated isotope shift is a quantity that 
remains unaffected by the A v coupling. 

As has been mentioned in earlier publications 0,0,0]) typical relativistic parametrizations 
predict a neutron skin in 208 Pb of about 0.3 fm. This should be contrasted with the case of 
non-relativistic forces of the Skyrme or Gogny type that yield values that are significantly 
smaller (between 0.1 and 0.2 fm). As the proton rms radius in 208 Pb is reproduced by both 
relativistic and non-relativistic models, the spread in the neutron rms radius in 208 Pb — owing 
to the different density dependence of the symmetry energy in different theoretical models — 
is about 0.1-0.2 fm. Within the relativistic mean-field models a spread (i.e., a theoretical 
uncertainty) of the neutron skin for a specific nucleus can be simulated in a controlled 
manner through the addition of the isoscalar-isovector coupling A v [Tol . I23I ] . Starting from 
the original NL3 parameter set and changing the isoscalar-isovector coupling constant from 
A v = (pure NL3 interaction) to A v = 0.025, the neutron skin in 208 Pb varies from 0.280 fm 
to 0.209 fm 0. Hence, in this formalism the induced theoretical uncertainty of the neutron 
radius of 208 Pb between these two extreme values of A v is approximately 0.07 fm. This spread 
for 208 Pb is only slightly larger than the expected experimental error in the measurement 
of the neutron radius in 208 Pb via the parity violating electron-scattering experiment at the 
Jefferson Laboratory (known as the PREX experiment) fl3j | . Thus we take the theoretical 
spread in the neutron radius of 208 Pb as a baseline to estimate a realistic uncertainty in 
those nuclei involved in APNC experiments. 

The neutron skins calculated with both of the relativistic parametrizations (namely, A v = 
and A v = 0.025) for stable and long-lived isotopes of Cs, Ba, Dy, and Yb are displayed in 
Fig. |H For a specific isotope, the shift between the two curves displayed in the figure gives 
the model spread (due to A v ) of the neutron skin of the given nucleus. For all the considered 
isotopic chains the neutron skins calculated with A v = (the stiffest symmetry energy) and 
A v = 0.025 (the softest symmetry energy) lie roughly along two parallel lines. In general, 
the effect of A v on the neutron skin is more prominent for the more asymmetric and heavier 
nuclei. This can be observed for the relatively long chain of Cs isotopes in Fig. 01 where 
the spread of the neutron skin varies approximately from 0.03 fm to 0.06 fm in passing from 



7 



123 Cs to 137 Cs. 

To the extent that in the RMF model used in this work the proton rms radius of a given 
nucleus is almost independent of A v , the theoretical spread of its neutron skin is therefore 
approximately equal to the one for its neutron rms radius. At present, the situation regarding 
our ignorance of the neutron skin of heavy nuclei is unlikely to improve because the existing 
experimental data cannot reduce the spread in the neutron radii of Cs, Ba, Dy, or Yb any 
further. Thus, we now examine whether the PREX experiment on the neutron radius of 
208 Pb may place important constraints on the neutron skin of APNC nuclei, including those 
deformed isotopes with an odd number of nucleons. Figure 03 evidences the strong linear 
correlation j^, 0, |23| between the neutron skin of 133 Cs, 138 Ba, 164 Dy, and 174 Yb and that 
of 208 Pb. The calculations include values for the isoscalar-isovector coupling constant A v 
ranging from to 0.025. Note that the correlation holds even for the case of deformed 
nuclei having an odd number of nucleons. With the culmination of the PREX experiment, 
the present theoretical spread in the neutron radius of 208 Pb of ~ 0.3 fm (relativistic and 
non-relativistic models altogether) 0, will be replaced by a genuine experimental 

error five times smaller; that is, Ar n ( 208 Pb) ~ 0.056 fm |23|. This experimental error in 
the neutron radius of 208 Pb will also be the experimental error in the neutron skin, as the 
proton radius is known with a much higher accuracy. If the expected experimental accuracy 
is attained, then the spread in the neutron skin of the several APNC isotopes will also be 
(indirectly) reduced. In turn, this result will impact favorably on APNC observables by 
appreciably reducing the nuclear structure uncertainty. 



III. RESULTS FOR APNC OBSERVABLES 

Parity violation in atomic systems arises from the interference between the parity con- 
serving electromagnetic interaction and the parity violating weak interaction. Although 
most of the binding energy of the atomic electrons comes from their attractive Coulomb 
interaction with the Z protons in the nucleus, the weak charge of the nucleus Qw induces 
a small correction in the binding energy and parity of the electronic wave functions. These 
small corrections can be measured experimentally. The experimentally measured quantity 
is the electric dipole amplitude between two electronic states that in the absence of parity 
violating effects would have the same parity, thereby making the amplitude vanish. The ob- 
servable related to the dipole amplitude can be parametrized as follows [13, E, HI IE IHj]: 



A P vc(N,Z) 



aZ)Qw(N, Z) = i{Z) Q S ™(N, Z) + AQ^ P (N, Z) 



(2) 



where £(Z) embodies the atomic structure contribution and Qw(N, Z) represents the weak 
charge of a nucleus of neutron number N and electric charge Z. 

The experimentally measured weak charge of the nucleus (without radiative correc- 
tions [20]) Qw(N, Z) differs from the standard model prediction, 



Q S ^(N,Z) 
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(3) 



by a nuclear correction factor that arises from the difference between proton and neutron 
one-body densities. That is, 
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q n ( P ) = / drf(r)p n (p)(r) 
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Note that in Eq. (J3J) # w is the Weinberg (or weak mixing) angle, while p n ^ in Eq. (JSJ) is 
the neutron (proton) density normalized to one, and f(r) is an electronic form factor that 
describes the spatial variation of the electronic axial- vector matrix element over the size of 
the nucleus. 

The atomic parity formalism "starts" by assuming that neutrons and protons have the 
same spatial distribution and then adds a nuclear-structure correction factor to quantify 
the differences between the actual neutron and proton densities. This is motivated by 
the fact that the proton density, which is needed to determine various parity conserving 
observables, is very accurately known. During the course of this presentation — and as has 
been done elsewhere (20, — an additional quantity that parametrizes the nuclear structure 
corrections will be extensively used. This quantity is defined as follows: 

Q^ cl (iV, Z) = -N(q n - 1) + Z(l - Asm 2 6 w )(q p - 1) . (6) 

Note that, in terms of this quantity, the weak charge of the nucleus may be written as 
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To compute the axial form facto r f( r) one expands the electronic Dirac wave functions 
in a power series about the origin [20| (the Coulomb potential is assumed to be that of 
a uniform nuclear charge distribution). Using this prescription yields the electronic form 
factor in closed form: 
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where R p is the cutoff radius of a sharp nuclear charge density. Effects from the finite nuclear 
size are computed from the ground-state expectation value of the time-like component of the 
weak vector current by considering the neutron and proton as point-like particles. Again, 
we assume sharp proton and neutron densities with rms radii equal to those predicted by 
self-consistent RMF models. Note that the relation between the sharp cutoff radii R n ( P ) and 
the corresponding rms radii r n ( p ) is given by 
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Estimates for the weak matrix elements q p and q n are then given by the following expressions: 
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Although in the present approximation q p does not depend on nuclear structure effects, the 
sensitivity of q n to nuclear structure uncertainties appears in the form of the rms proton 
radius — which is well known — and the poorly known neutron skin t. 



A. Prospects for isotope ratios analyses 



The measurement of the weak charge in APNC experiments is plagued by theoretical 
uncertainties in both atomic and nuclear structure. Fortunately, uncertainties in atomic 
structure may be eliminated, or at least considerably reduced, by stud ying parity violation 
along a chain of isotopes and taking ratios of APNC measurements jl4L flsL Eol . 21, 2^ . 
According to Eq. (J2J), the dependence of the parity violating amplitude Ap^c on the atomic 
theory contribution £(Z) will cancel out in the ratio of two measurements performed in two 
different isotopes of the same element, provided £,(Z) does not change appreciably along 
that isotopic chain. Thus, we now address the role of the remaining — nuclear structure — 
uncertainty in APNC observables along the isotopic chains of Cs, Ba, Dy, and Yb. The 
observables of interest are expressed in the form of the following two ratios: 



Q W (N',Z)-Q W (N,Z) 
Q W (N',Z) + Q W (N J Z) 
Q™(N>,Z) + AQ n w p (N',Z) 



Q S ^(N,Z) 



AQ^ P (N, Z) 



Q s /(iV', Z) + AQ n w p (N>, Z) + Q s /(iV, Z) + AQ n w p (N : Z) 



Q™(N',Z) + Qff*{N',Z) 



nucl / 



(N,Z) 



/Onucl 



(N,Z) 



QSM(iV', Z) + Q^ C \N' } Z) + QSM(AT Z) + Q^(N, Z) 



(12) 



and 



tu 



Qw(N',Z) 
Qw(N,Z) 



Q™(N',Z) + AQ n w p (N>,Z) 
Q™(N,Z) + AQ n w p (N,Z) 
Q™(N>,Z) + Q™ C \N',Z) 
Q S ™(N,Z) + Q™ C \N,Z) ■ 



(13) 



The approximate sign («) in the above two equations follows from the assumption that the 
overlap q p appearing in Eq. (J2J) remains constant along the whole isotopic chain. As can be 
realized from Tables HT1 and |llll (see below), this is an excellent approximation. 

It has been argued in Ref. [22j that corrections to standard model predictions induced by 
new physics or uncertainties in nuclear structure are essentially the same whether one uses 
1Z\ and 1Z2- Thus, we focus here exclusively on TZ\. Using Eqs. Q and (0) we can write the 
ratio TZ\ approximately as 



Tlx 



Q S ^(N',Z) 



Q s ^(N,z) 



Q™(N>,Z) + Q™(N,Z) 



1 + 



N' 
A~N 



[q n (N',Z)-q n (N,Z)} 



(14) 



where N' (N) is the largest (smallest) neutron number and AN = N' — N represents the 
difference in neutron number between the two extremes of the isotopic chain. Further, 
Ref. 0] also established that a significant test of the standard model on the basis of isotope 
ratios requires a measurement of IZi with an accuracy better than 0.1%. Otherwise, a 
less precise determination of isotope ratios would not be worthwhile to compete with the 
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sensitivity to new physics of the presently known data for the cesium isotope 133 Cs. The 
measurement of the electron-nucleon parity violating effect in 133 Cs to 0.35% accuracy [3^ 
remains to date as the paradigm of APNC experimental precision. 

Because atomic uncertainties have been eliminated from the ratio 1)140. the remaining 
uncertainties in TZi reflect the known accuracy in the neutron and proton rms radii. While 
proton densities have been determined from electron-scattering experiments with remarkable 
accuracy, knowledge of the neutron densities to a comparable level of accuracy is lacking. 
Thus, the main nuclear structure uncertainty in the isotopic ratio TZi comes from our poor 
knowledge of the neutron radii (or neutron skin) of heavy nuclei. Indeed, by inserting 
Eq. 1)11)1 into Eq. 1)140. the relative uncertainty in IZi may be approximated by 



Tlx 525 V ; AN 



±(N',Z)--(N,Z) 



(15) 



In Sec. |H] we showed that a spread in the neutron radius of 208 Pb is obtained by using 
parametrizations with isoscalar-isovector couplings A v = and A v = 0.025, which is represen- 
tative of the experimental error (~ 1%) expected from the Jefferson Laboratory experiment. 
Note that the quoted error in the neutron radius of 208 Pb is also the error in its neutron skin 
as the proton radius is known with a much higher accuracy. Inasmuch as a linear correlation 
between the neutron skin of 208 Pb and that of the APNC nuclei has been established (see 
Fig. |SJ), an estimate of the theoretical uncertainty in the weak charge can be obtained from 
the spread in the neutron skin of the APNC nuclei using the A v = and the A v = 0.025 
parameter sets. 

Let us outline the procedure we shall follow to determine the relative error in TZi (namely, 
SIZi/TZi). First, we calculate the neutron skin t = r n — r p in a certain isotopic chain using 
the A v = parametrization. Since r p is roughly similar for all the isotopes of the chain, the 
fractional neutron skin variation along the chain is estimated as 

t(N', Z) - t(N, Z) 
(At) Av=0 = r~, . , (16) 



(r P ) 



A v =0 



where (N', Z) refers to the heaviest member and (N, Z) to the lightest member of the isotopic 
chain under consideration, and (r p ) is the average proton rms radius of the nuclei in the 
chain. Next, we compute the corresponding fractional neutron skin variation (AF)a v= o.o25 
for the softer A v = 0.025 parameter set. Finally, the model spread 2^| is obtained as the 



difference between the fractional neutron skin variations calculated in the two parameter 
sets. That is, 

S(At) = (Ai) Av=0 - (At) Av=0 . 25 • (17) 

In Table H] we display the relative uncertainty 51Zi/TZi in the ratio TZi for each one of 
the isotope chains (Cs, Ba, Dy, and Yb) together with the different components required 
to compute it. These are the largest neutron number of the isotopic chain N', the largest 
difference in the neutron number of the isotopic chain AN, the average proton rms radius 
(r p ), and the model spread 6 (At). The latter turns out to be a small quantity, not only 
in the present relativistic calculations but also in calculations with non-relativistic nuclear 
energy functionals (see Ref. [i^l and references quoted therein). The reason resides in the 
approximately linear relationship of the neutron skin along a chain of isotopes with the 
neutron number. While the model dependence in the value of the intercept is obvious (see 
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Fig. HJ), the slope is only mildly model dependent. And it is the difference in slope that the 
model spread 6 (At) is sensitive to. 

It is seen that the estimated relative uncertainty of the APNC ratios TZ\ in the four 
isotopic chains considered is around 0.25%-0.35%, with a slight increase with increasing 
atomic number. The quoted uncertainty appears somewhat larger than the 0.1% value that 
would be desirable to extract new physics limits beyond the present reach of the 133 Cs 
experiment. Note that the neutron skin along a given isotopic chain varies (almost) linearly 
with neutron number (see Fig. 0J). This makes the model spread 5 (At) roughly proportional 
to the neutron difference AN. As the combination 5 (At) / AN is therefore nearly independent 
of A7V" |38j], it is unlikely that the precision in TZi can be improved by enlarging the range 
of the neutron difference AN between the two extremes of the isotopic chain [see Eq. ([15)1] . 

Table |U also displays the model spread [<5(Af)] ai % that would be required to reach the 
sought-after relative uncertainty of 0.1% in the isotopic ratio Hi. In all cases this number is 
smaller than the value provided by the RMF model. Recall that we have tuned the effective 
interaction by means of A v to mimic the purported 1% accuracy of the PREX experiment. 
Hence, it appears that under the present — and future — situation, nuclear structure uncer- 
tainties affecting the variation of the neutron distribution along an isotopic chain are too 
large to make 7Zi a useful probe for physics beyond the standard model. However, note that 
the differences are not dramatic: a factor 2.5 for Cs and Ba, 2.75 for Dy, and 3 for Yb. Thus, 
second generation experiments and/or novel facilities may significantly aid in this quest. 



B. Prospects for single-isotope analyses 

For most of the APNC nuclei — particularly Cs and Fr, with one and no stable isotopes, 
respectively — one cannot envisage obtaining precision data on several members of the iso- 
topic chain from experiments at the present time. For these cases, and to provide useful hints 
for future experiments, accurate estimates of the uncertainty in the neutron skin are essen- 
tial. Thus, we now compute the weak charge — including nuclear and one-loop electroweak 
radiative corrections — for some Cs isotopes using both relativistic (A v = and A v = 0.025) 
parametrizations. The relevant formulas to be employed are given by Eqs. (jSJ)-© suitably 
modified by the inclusion of radiative corrections derived by Marciano and Rosner j40j; that 
is, 



Q^(N, Z) = (0.9857 ± 0.0004) 
Qw c \N, Z) = (0.9857 ± 0.0004) 



N + Zll- (4.012 ±0.010)x 



(18) 



N(q n -1) + Z[1- (4.012 ±0.010)x)(g p - 1) , (19) 



where we have adopted the shorthand notation 



x = sin 6 



w 



0.2323 



(20) 



Notice that corrections due to new 
these may be found in Refs. lil 



hvsics have not been included in the above expressions; 

To be able to compare more precisely our results with those in the available literature, 
the weak matrix elements q p and q n appearing in Eq. (JT§j) have been slightly modified from 
those in Eqs. (jlOj) and (jllj) . We now take into account the corrections due to the finite size 
of the proton and neutron charge distributions, which originate additive correction terms q c v 
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and q c n to Eqs. (fTUjl and (fTTJl . respectively. These terms are given by [2 



5Vi(p) 



Jdr 1 - <r 2 )^ (p) /(r) V^/Q^), (21) 



where (t" 2 )/^,) are the intrinsic nucleon weak rms radii, p c n ^ are the density distributions 
of nucleon centers, and Q^ p ^ are the nucleon weak charges. Assuming uniform nucleon 
distributions, we can cast the corrective terms q°^ due to the finite size of the nucleon 
as 

<--™'t?("-^). 

* = -< z ^( - M -£jf)- (23) 

where i? denotes the sharp cutoff charge radius (which now includes the finite size of the 
proton) and a contribution of 0.14 fm 2 from the strangeness radius of the nucleon is included. 

The analysis of the data of a clean measurement by Wood et al. [33] of the amplitude of 
the parity non-conserving transition between the 6s and 7s states of 133 Cs, the only naturally 
occurring cesium isotope, established a value of Q^ P = — 72.06(28) e xp(34) t heor for the weak 
charge of this element [15]. This experimentally extracted result differed by 2.5 o from the 
standard model prediction with radiative corrections of = —73.20 ± 0.13 [Ref. 40] and 
Eq. (J18)) ]. Let us mention that a is evaluated by adding in quadrature the experimental 
and theoretical errors quoted for Q^f. The excitement over this apparent disagreement 
has faded as a consequence of recent reports that demonstrate that the deviation between 
the experimental results and the standard model prediction can be removed bytaking into 
account self-energy and vertex QED radiative corrections j41, 42| (see also Ref. for an up- 



to-date review of the status of measurements and calculations of parity violation in atoms). 
The newly reported values of the weak charge for 133 Cs — which reconcile the experiment 
with the standard model prediction — are Q^f = — 72.71(29) cxp (39) t heor |4jJ and = 



— 72.90(28) ex p(35)theor |42|. It should be noted that in the extraction of the experimental 
values of the weak charge Q^ p , one assumes equal proton and neutron densities (normalized 
to unity), so that the overall nuclear correction (namely q p ) factors out. An efficient way to 
account for the difference between neutron and proton densities — and thus between q n and 
q p — has been indicated in Eq. (j3j). Including radiative corrections, this nuclear structure 
correction factor becomes 

AQ^ P = 0.9857iV (\ - ^ , (24) 

and the total weak charge of the nucleus is computed as Qw = Qw 1 + AQ^7 P . 

Results for the weak charge together with other relevant quantities, such as q n ^ [including 
the intrinsic structure corrections (|22|) and (J23|) ]. the proton rms radius r p , and the ratio r n /r p 
of neutron to proton rms radii are displayed in Table HP for Cs isotopes. As expected, we find 
that q p is constant, model independent, and within 0.04% of the sharp-radius value expressed 
in Eq. (jlOjl . However, both the changes and the model dependence of the neutron radii along 
the isotopic chain are noteworthy. Clearly, the ratio r n /r p increases with increasing neutron 
number N. Further, for a given neutron number N, r n /r p decreases with increasing A v , 
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owing to the softening of the symmetry energy. Hence, q n decreases with increasing neutron 
number and with a decreasing A v . This ultimately leads to an increase in the nuclear part 
of the weak charge Q^ d [see Eq. ()19|)]. In summary, while q p is largely independent of A v , 
q n shows a moderate increase with A v . Thus, through variations in A v one has a margin to 
adjust the nuclear contribution to the weak charge. In particular, with the proposed model 
dependence a change of 0.062 is induced in Q^} cl for the 133 Cs atom. 

The values of the nuclear structure correction factor AQ^ P for the studied Cs isotopes, 
calculated from (j24j) . have been tabulated in Table ITT1 For the specific case of 133 Cs, we find 
a nuclear structure correction factor of AQ^ P = 0.342 (0.278) for the A v = (A v = 0.025) 
parameter set. Addition of this nuclear correction term to the standard model prediction 
with radiative corrections (—73.197, third column of Table ITT]) ends in a total weak charge 

ft^Cs) M2.855, ifA v =0, 

WV ; [-72.919, ifA v = 0.025, V ; 

which is in satisfactory agreement with the measured values according to the recent revi- 
sions H Hi . 



It has been discussed in the literature [24], |25j that AQ^ calculated with RMF 
parametrizations is almost twice as large as the corresponding values predicted by Skyrme 
interactions 21]. This is a reflection of the stiffer symmetry energy of the relativistic models. 
For example, by using the sharp-cutoff-radius approximation of Eqs. (fTU|) and (|TT|) one can 
write 

AQ n w p = 0.9857N { -^^-. (26) 

Because the above nuclear correction factor is directly proportional to the neutron skin 
t, RMF parametrizations with stiffer symmetry energies — and consequently larger neutron 
skins — yield larger values for AQ^ P relative to Skyrme interactions. Equation (|26|) may 
also be used to estimate the model spread in AQ^ P for 133 Cs. This is given by 

5 (AQ n w p ) = 0.9857AT ^£ ^ 6 ( -) « 0.064 , (27) 



q p 525 \r p/ 

where the model spread for t/r p has been derived from the r n /r p values calculated with 
A v = and A v = 0.025 for 133 Cs (listed in Table ITT]) . A straightforward estimate of the spread 
computed as the difference of the values of AQ^y P between the A v = and the A v = 0.025 
parameter sets yields 

5 {AQw P ) « 0.342 - 0.278 = 0.064 , (28) 

which compares well with our previous value. In Ref. |43[ it has been argued that for 
the purpose of significantly reducing the nuclear structure uncertainties, relative to those 
arising from radiative corrections, the neutron radius of 133 Cs must be known within 2%. The 
promise of a 1% measurement of the neutron radius of 208 Pb at the Jefferson Laboratory [l3| . 
together with theoretical correlations of the type displayed in Fig. will be of considerable 
help in reaching the desired goal. 

C. A glimpse into the future: Fr isotopes 

Another strategy proposed as a means for attempting to improve the accuracy in APNC 
experiments is the utilization of high-Z atoms that enhance parity violating effects owing 
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to both the increase in the number of nucleons and the increased electronic density in the 
neighborhood the nucleus. Further, alkali atoms have the simplest atomic structure, thereby 
minimizing theoretical uncertainties. This makes francium, the heaviest "simple" atom, an 
attractive candidate for APNC experiments. Indeed, the increased number of nucleons and 
the increased electronic density in the re gion of the nucleus makes the APNC effect in Fr 
some 15 to 20 times larger than in Cs Eij ]. Unfortunately, francium has no stable 
naturally occurring isoto pes . Although remarkable experimental progress has already been 
achieved in the field [i^, 0, |46| , the production of intense beams of radioactive Fr isotopes 
represents an exciting challenge for the future of the rare-isotope facilities all over the world. 

We present estimates for APNC observables in francium along an isotopic chain that 
ranges from 207 Fr to 225 Fr and that includes the two isotopes with the longest half-life (of 
about 20 min); these are 212 Fr and 223 Fr. In Fig. [Hjwe display the neutron skin, charge radii, 
binding energy per nucleon, and deformation for six isotopes in the chain computed with 
the two extremes of the values of the isoscalar-isovector coupling constant used in this work 
(A v = and A v = 0.025). As established earlier for 174 Yb in Fig. Q the calculated values 
for the charge radii, binding energies, and deformation show little model dependence and 
reproduce the existing experimental data (where available) 0, 0| rather well. As expected, 
it is only the (experimentally unknown) neutron skins that are sensitive to the value of A v . 
Note that with the exception of the semi-magic nucleus 213 Fr, the ground state of all other 
isotopes is found to be deformed. 

As was done for the Cs isotopes in Table |HJ we display in Table 11111 the calculated weak 
charge for the six considered Fr isotopes, alongside other relevant observables. According 
to Eq. (|18|) . the standard model predictions for the weak charges of 212 Fr and 223 Fr result 
in QSM = _ 117 . 380 ± 0.20 and -128.223 ± 0.20, respectively, where the uncertainty in 
has been estimated in a similar way to that of 133 Cs in Ref. |4£j. The nuclear corrections 
to the weak charge of 212 Fr ( 223 Fr) calculated from q n and q p through Eq. (|2*3j) are given as 
follows: 



These results suggest that nuclear-structure uncertainties in our model, arising from dif- 
ferences in the shape of the neutron density relative to that of the proton, are 0.266 and 
0.347 for the francium isotopes 212 Fr and 223 Fr, respectively. The numbers are similar to the 
uncertainties 0.264 and 0.345 that one would obtain with the sharp-radius approximation 
of Eqs. and 

Thus, with the expected precision to be attained by the PREX experiment, which we 
mimic here through a change in A v , our predictions for the weak charge Qw = Qw 1 + ^Q\v P 
of the 212 Fr and 223 Fr isotopes vary between Q w = -116.209 (A v = 0) and -116.475 
(A v = 0.025) for 212 Fr, and between Q w = -126.639 and -126.986 for 223 Fr. Note that the 
uncertainty due to A v should be augmented by a ±0.20 error from Qyy~. 

IV. CONCLUSIONS 

Motivated by the prospects of a high-precision (1%) measurement of the neutron radius 
of 208 Pb at the Thomas Jefferson Laboratory, we have examined the impact of such experi- 
ment on atomic parity non- conserving observables. While such a measurement will have far 
reaching consequences in fields as diverse as nuclear structure, flavor violations in the strong 




(29) 
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interactions, and nuclear astrophysics, we have focused in this contribution exclusively on 
APNC observables. 

Our theoretical framework is based on the highly successful relativistic NL3 parameter 
set, suitably modified by the addition of an isoscalar-isovector term (A v ). The virtue of such 
a term is that it enables one to change the neutron radius of heavy nuclei — which at present 
is poorly known — without compromising the success of the model in reproducing a variety 
of ground state properties of nuclei throughout the periodic table. Relative to an earlier 
study based closely on this formalism, the present study improves on it in one essential 
aspect: the inclusion of deformation and pairing correlations. Without adding these effects, 
the predictions for the structure of most nuclei along the isotopic chains of relevance to the 
APNC program would be unreliable. As a result, in the present work we have mapped the 
neutron radii of the relevant isotopic chains as a function of both N and A v . As established 
by others before us, we have found a strong correlation between the neutron radius of 208 Pb 
and the neutron radius of APNC nuclei — even for the case of odd-even and odd-odd nuclei. 
Employing this correlation, we have used a range of values for A v that closely matches 
the 1% (or ~ 0.056 fm) uncertainty expected to be achieved by the Jefferson Laboratory 
experiment. 

With this information at hand, we proceeded to study the impact of such a 1% measure- 
ment on two different combinations of APNC observables. The first set of observables (Hi 
and 7^2) are formed from ratios of weak nuclear charges along isotopic chains. The merit 
of such observables is that the ratios are largely insensitive to uncertainties in atomic struc- 
ture, leaving nuclear uncertainties — in the form of the neutron skin — as the main source of 
theoretical error. The second set of observables involves a direct determination of the weak 
charge of various alkali metals, such as cesium and francium. Unfortunately, these elements 
have very few stable or long-lived isotopes, so in their case the accurate determination of 
ratios of weak nuclear charges is more difficult. 

In the case of the isotopic ratio Hi, it has been claimed that a significant test of the 
standard model requires a determination of Hi to better than 0.1%. This precision would 
be required for 1Zi to supersede the 133 Cs experiment as the most stringent test for new 
physics within the APNC program. Unfortunately, our results indicate that the projected 
1% accuracy in the measurement of the neutron radius of 208 Pb at the Jefferson Laboratory 
appears unlikely to translate into the required 0.1% (or lower) uncertainty in Hi. Instead, 
we have established an uncertainty in Hi that is two to three times larger (of the order 
of 0.25-0.35%). Although the Jefferson Laboratory experiment is unlikely to achieve the 
desired accuracy, it is plausible that second-generation experiments may reach this goal. 

In the case of the weak nuclear charge of the alkali metals, where few stable or long- 
lived isotopes exist, one must consider theoretical uncertainties arising from both atomic 
and nuclear structure. Although it is no longer possible to fully eliminate uncertainties 
in atomic structure (as was done for Hi) alkali metals at least enjoy the simplest atomic 
structure. As far as the nuclear structure is concerned, the uncertainty is fully subsumed 
into a single factor: the neutron skin of the nucleus of interest. It has been argued that 
to reduce significantly the nuclear-structure uncertainties relative to radiative corrections, 
the neutron radius of 133 Cs must be known with an accuracy of 2% or better. Similar 
accuracy should be expected for the case of other alkali metals. Indeed, in anticipation of 
the commissioning of rare-isotope accelerators all over the world and the ongoing advances in 
the field of production and magneto-optical trapping of radioactive atoms, we have computed 
the weak nuclear charge of the two longest lived francium isotopes: 212 Fr (20 min) and 223 Fr 
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(21.8 min). Our assumed uncertainty in the neutron radii of 212 Fr and 223 Fr translated 
into a theoretical uncertainty in the value of their weak nuclear charges of 0.2% and 0.3%, 
respectively. 

In summary, we have studied the impact of a high-precision measurement of the neutron 
radius of 208 Pb on atomic parity non-conservation experiments. However, the relevance of 
such a measurement on a plethora of other research areas — such as nuclear structure, flavor 
violations in QCD, and neutron-star structure — has been strongly emphasized here and 
elsewhere. Thus, the neutron radius of 208 Pb stands as one of the most fundamental nuclear- 
physics quantities yet to be accurately determined. This unfortunate situation should be 
promptly corrected. 
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TABLE I: The relative uncertainty 51Z\/1Zi in the APNC isotope ratio 1Z\ and the various com- 
ponents needed to evaluate it, as described in the text. The last row denotes the model spread 
5 (At), Eq. (|17j). that would be required to achieve 5TZ\/TZi = 0.1%. 



Observable 


Cs 


Ba 


Dy 


Yb 


Z 


55 


56 


66 


70 


N' 


82 


82 


98 


106 


AN 


14 


8 


8 


8 


(r P ) 


4.752 


4.766 


5.134 


5.254 


5(At) 


0.0059 


0.0033 


0.0022 


0.0022 




0.0025 


0.0025 


0.0027 


0.0035 


[*(At)] .i% 


0.0024 


0.0013 


0.0008 


0.0007 



TABLE II: A variety of quantities of relevance to APNC (as defined in the text) for four Cs isotopes 
(Z = 55), including the sole stable element 133 Cs. The experimental value for the weak char ge o f 
this isotope is given by Q™ p = -72.71(29) ex p(39) theor Q and Q^ p = -72.90(28) exp (35) theor |3. 



A v 


N 




/onucl 




Q P 




r p 


AQ n w p 


Qw 


0.000 


76 


-71.226 


3.284 


0.95410 


0.95803 


1.04103 


4.74141 


0.308 


-70.918 




78 


-73.197 


3.400 


0.95377 


0.95803 


1.04568 


4.74551 


0.342 


-72.855 




80 


-75.169 


3.512 


0.95350 


0.95803 


1.04949 


4.74979 


0.373 


-74.796 




82 


-77.140 


3.622 


0.95327 


0.95804 


1.05277 


4.75903 


0.402 


-76.738 


0.025 


76 


-71.226 


3.228 


0.95484 


0.95803 


1.03062 


4.75244 


0.249 


-70.977 




78 


-73.197 


3.338 


0.95457 


0.95804 


1.03449 


4.75720 


0.278 


-72.919 




80 


-75.169 


3.444 


0.95436 


0.95804 


1.03742 


4.76113 


0.302 


-74.867 




82 


-77.140 


3.549 


0.95418 


0.95804 


1.04007 


4.76972 


0.325 


-76.815 
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TABLE III: A variety of quantities of relevance to APNC (as denned in the text) for six Fr isotopes 
(Z = 87). Francium has no stable isotopes. 



A v 


N 






Qn 


q P 


r n /r p 


r p 




Qw 


0.000 


120 


-112.451 


12.650 


0.88788 


0.89510 


1.02993 


5.51834 


0.954 


-111.497 




125 


-117.380 


13.361 


0.88659 


0.89510 


1.03724 


5.53761 


1.171 


-116.209 




126 


-118.366 


13.502 


0.88636 


0.89510 


1.03857 


5.53974 


1.213 


-117.153 




132 


-124.280 


14.294 


0.88544 


0.89511 


1.04402 


5.60901 


1.406 


-122.874 




136 


-128.223 


14.866 


0.88454 


0.89511 


1.04921 


5.64364 


1.584 


-126.639 




138 


-130.194 


15.171 


0.88397 


0.89512 


1.05246 


5.66462 


1.693 


-128.501 


0.025 


120 


-112.451 


12.453 


0.88955 


0.89510 


1.02062 


5.52739 


0.734 


-111.717 




125 


-117.380 


13.123 


0.88853 


0.89510 


1.02641 


5.54760 


0.905 


-116.475 




126 


-118.366 


13.250 


0.88839 


0.89510 


1.02723 


5.55034 


0.932 


-117.434 




132 


-124.280 


14.015 


0.88759 


0.89511 


1.03199 


5.61457 


1.094 


-123.186 




136 


-128.223 


14.556 


0.88686 


0.89511 


1.03624 


5.65014 


1.237 


-126.986 




138 


-130.194 


14.844 


0.88637 


0.89512 


1.03904 


5.67325 


1.329 


-128.865 




0.005 0.01 0.015 0.02 0.025 



A 

V 

FIG. 1: The dependence of the ground-state charge radius, binding energy, quadrupole deforma- 
tion, and neutron skin on the isoscalar-isovector coupling A v is exemplified for the 174 Yb nucleus. 
The observables have been normalized to the values of the NL3 interaction (A v = 0). 
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FIG. 2: Charge radii calculated with the models considered in the text are compared with the avail- 
able experimental data 0] in isotope chains of interest for APNC studies. The actual experimental 
error bar for 156 Dy and 158 Dy has been reduced by a factor of 10 for display. 
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FIG. 3: Isotope shifts of the charge radii of cesium and francium nuclei. The experimental data 
are from Ref. 



21 



0.25 - 



0.20 - 



0.15 - 



0.10 - 



168 170 172 174 176 156 158 160 162 164 



S-H 
I 



0.25 



0.20 - 



0.15 - 



0.10 - 




•— 


— • 


A = 

V 







— ♦ 


A = 

V 


0.025 





FIG. 4: Variation of the neutron skin in four isotope chains of possible relevance to APNC isotope 
ratios measurements, for the two effective field theory models considered in this work. Note that the 
change in the neutron skin of 208 Pb predicted by these two models, otherwise accurately calibrated 
for binding energies and charge radii, mocks up the uncertainty that will be left after the projected 
precision measurement of the neutron radius of 208 Pb at the Jefferson Laboratory 
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FIG. 6: Calculated ground-state properties (neutron skin, charge radius, binding energy per 
particle, and quadrupole deformation) of some francium isotopes. Experimental data are shown 
for comparison when available 0, |4tJ ■ 
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